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Entanglement [l| has been recognized as the essential 
element of quantum information science. This is due to 
the fact that the nonlocal quantum coherent nature of en- 
tanglement can be used as a resource for implementing 
quantum information protocols "^l. Recently, the concept 
of entanglement has also been introduced to the inves- 
tigation of quantum phase transition(QPT)[^ which is 
induced by quantum fluctuations and therefore can oc- 
cur even at zero temperature. More precisely, QPT can 
be identified as the appearance of non-analyticity in the 
ground state energy. For the case of spin-boson model 
treated in this work, a phenomenon closes to QPT known 
as Ground State Instabilities (GSI) also occurs in finite 
system. One will see that non-analyticity arises due to 
level crossing^ which indicates the instability of ground 
state. Similar problems of QPT in the weak coupling 
and thermodynamical limit have been studied by many 
authors ,5|. In addition, the related problem of entan- 
glement in the so-called Dicke model(DM)'6'] has also at- 
tracted much attention|2| recently. One of the interesting 
results of the Dicke model is the sequence of GSI in arbi- 
trary finite-atom systemj^ which has been overlooked in 
the thermodynamical limit. More interestingly, at these 
infinite sequential instabilities, Buzek et al. show that 
there are corresponding discontinuities appearing in the 
ground-state entanglement of the reduced atomic system. 
However, it has been pointed out by K. Rzazewski and K. 
Wodkiewicz^J that in ^ gauge invariance is spoiled by 
not including the term of the minimal coupling hamil- 
tonian. Furthermore, they also pointed out that without 
the contribution the hamiltonian is unbounded from 
below as the coupling goes to infinity. It is obvious that 
any two-level atomic system is isomorphic to a spin-1/2 
system. Therefore the Dicke model can be identified as 
a spin-boson interacting system. Certainly for the spin- 
boson system, there is no requirement of gauge invari- 
ance. Furthermore, by keeping finite coupling it seems 
that the unbounded problem can be avoided. However, 
by requiring finite coupling, one can only allow the in- 



vestigation of finite number of ground state transitions 
instead of the infinite transitions in Even though 
with such restriction, the correlation between entangle- 
ment and GIS can still be addressed. In order to un- 
derstand the relation between ground state instability 
and concurrence, exact analytical results are needed for 
gaining insight. Here, we discuss the generic spin-boson 
model by introducing a parameter r which is the detun- 
ing of the boson mode frequency. Moreover, for more 
realistic consideration, we also study the ground state 
instabilities in the multi-mode model. In this work, we 
rigorously show the existence of sequential ground state 
instabilities for arbitrary detuning in the two-spin sys- 
tem. In contrast to the resonant case where GIS and 
concurrence are strongly correlated, such detuning effect 
leads to the disconnection of ground state instability and 
concurrence. By including extra modes and tuning the 
frequency, it is found that the ground state can become 
entangled even in the weak coupling regime and the en- 
tanglement is enhanced comparing with the mono-mode 
on-resonant spin-boson model(SBM, from now on, the 
term SBM denotes the mono-mode on- resonant model). 
The plan of the paper starts by introducing the single- 
mode spin-boson model, and the exact spectrum is then 
presented for two-spin system. In the next section we 
will show that ground state instability is a generic phe- 
nomenon of spin-boson model. Section HI provides the 
analysis of ground state entanglement by calculating the 
concurrence. In this section, we establish the fact that 
ground state instability is not directly correlated with the 
analyticity of concurrence. The study of ground state in- 
stabilities of the two-mode model is presented in section 
IV. By including extra modes, for two spins, it is found 
that there exists a region of detuning where enhanced 
ground state entanglement can be obtained even in the 
weak coupling regime. Finally, a brief summary is given 
in the last section. 



I. THE SINGLE-MODE MODEL AND ITS 
SPECTRUM 
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To begin with, we discuss the general method to solve 
the iV-spin model. The system is N spins interacting 
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with a mono-mode boson field. The Hamiltonian of the 
total system in the interaction hamiltonian is given by 

H — ojqJz + t^a^CL + gJ+a + g*J-a^ (1) 

where = ^ "^jLi '^fj ol can either be — } for raising 
and lowering operations or {a;,?/,z}, 0(0^) is the boson 
annihilation(creation) operator, ct" are the Pauli ma- 
trices. [jJo is the level spacing of the spin and w indi- 
cates the frequency of the boson mode. We have assumed 
that these spins couple to the boson mode with the same 
strength g. This detuned spin-boson model(DSBM) can 
in principle be solved ex actly JUl. In this work, we ex- 
tend the method of Swain jllj to diagonahze DSBM. The 
Hamiltonian can be separated by = Hq + Hj: 

Ho = J, + a'^a (2) 
Hj = ra) a + njj^a + K* J^a) (3) 

where r = uj/ujq — 1 and k = g/ujo. The parameter 
r(— 1 < r < 00) is related to the detuning which is usu- 
ally defined by w — wq in quantum optics. Hq is the 
so-called excitation operator Q. To obtain the spectrum, 
one uses the fact that {H, Hq, Hj, J^} form a maximally 
compatible set, where J"^ = + Jy + . For a iV-spin 
system, we focus on j = N/2 which is relevant to the 
ground state discussions. Due to the commutative J^, 
the matrix of H is automatically block diagonal by each 
j in the basis of Hq which is denoted by {\j,iTi)A\'n)p}. 
\j,m)A and \n)p are the spin states and photon number 
states respectively. {A = m + n} are the eigenvalues of 
{IjiTTT^Aln^p}- It is noted that by excluding ra^a our 
Hq = Jz + a'' a is a. parameter free operator and so does 
its eigenvalues. This approach helps to ease the count- 
ing of degeneracy of Hq. Hence, for fixed A, there exists 
degenerate subspace such that the diagonalization of H 
reduces to diagonalize finite matrices of Hj. The eigen- 
states of H are denoted by |j. A, h). Then 

H\j,X,h) = Exh\jA,h) (4a) 
Hi\j,X,h) = h\j,\,h) (4b) 

\j,X,h) = ^ A^■''''^\J,\^^)A\^)p (4c) 

Exh ^ \ + h. (4d) 

Note that, with arbitrary detuning, the energy eigenval- 
ues depend not only on k, but also on r and one should 
expect some new results due to these parameters depen- 
dence. The detail form of Hj for arbitrary N are given 
in Appendix A. Since the dimension of Hj becomes big- 
ger as the excitation number increases, most spectrum 
can only be obtained numerically for iV > 3. However, 
for two-spin, the full spectrum with j — I can be ob- 
tained(we further neglect the index h since only the eigen- 
value of Hj which is a decreasing function of k is needed 



for fixed A.): 

El = -1 (5a) 
Eq = i{r- v/8kP + r2} (5b) 

Ex = A + Ar-^V3^cos{|-^} (5c) 

ax = {4X + 2)\k\^ + r"^ 

_i ,3\/3K^r 
ipx = cos { 

where A runs from 1 to infinity. 



II. SEQUENTIAL GROUND STATE 
INSTABILITIES IN DSBM 

Usually, the eigenenergies of a quantum system are an- 
alytic functions of the coupling constant k. However, 
there is a possibility that when H{k) = Hq + kHj and 
[Ho,Hi] = such that Hq and Hj can be simultane- 
ously diagonalized and therefore the eingenfunctions are 
independent of k even though the eigenvalues vary lin- 
early with Kf4^. As a result, when one of the excited 
state is crossing with the ground state at some criti- 
cal value K = k, non-analyticity appears in the ground 
state energy. Such level-crossing phenomenon is called 
ground state instability(GSI) which also happens in the 
system considered in this work. The level crossing of 
SBM can be illustrated easily from the eigenstates with 
and 1 excitations. From Eq.(A3), it is clear that as 

K < + r)/N the energy En \s less than En so 

2 ^ 2 +-■- 

the ground state is the one with zero excitation. However, 

for K > + r)/N, the ground state is replaced by the 

one with 1 excitation since En_,-, < En_- At the critical 

value Kn_ = \/(l -f r)/N, the excitation number changes 

2 

discontinuously from A = — ^toA^— y + Obvi- 
ously, when r — 0, Kn_ reduces to the on-resonant result 

2 

in Ref. 8]. Buzek et al have shown numerically that the 
ground state energy is non-analytic and the level crossing 
occurs in sequence: {Ei Eq ^ Ei ^ E2 ^ ■ ■ ■}■ Due 
to the fact that the spectrum for two-spin can be obtained 
in closed form, we can provide an analytic proof for these 
sequential ground state transitions for all r if the follow- 
ing conditions are satisfied: (av>0,— l<r<c» and 
A > —1): (i) {Ex} are monotonic decreasing functions, 
except for A = — 1. (ii) For all A, /(k, r, A) = Ex+i — Ex, 
f{K, r, A) is a monotonic decreasing function with oppo- 
site signs at small and large k. With kx denoted the value 
of level crossing which is determined by Ex+i = Ex, 
we have {in) {Ex+2 > Ex+i}\K=iix- Due to the abso- 
lute square of k in Eq.(5), one may choose k > with- 
out losing generality. The first condition guarantees the 
eigenenergies of different A involved in the ground state 
level crossing at different coupling strength. The second 
condition ensures that there is only one crossing between 
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Ex and Ex+i- If Ex+2\ii^ is larger than Ex+ila^, then 
the crossing kx+i which is determined by the equation 
Ex+i = Ex+2 raust be larger than kx- Therefore, these 
three conditions together ensure GSI occur in sequence. 
The detail proof is given in Appendix B. 

It is interesting to point out that by adjusting the de- 
tuning parameter r, it is possible to have GSI in the small 
coupling regime. We recall the fact that, for on-resonant 
case(r = 0), the first ground state transition occurs at 
Kj = \f\j2. However, with negative r(u) < loq), the crit- 

'1/2 resulting with 



ical Ki = •\/(l -f r)/2 is smaller than 
GSI at weak coupling. On the contrary, if lo is larger than 
ujQ, then one needs a strong coupling strength to obtain 
GSI. Let 

'^f J '^i ^n-'i '*f the ith GSI critical couplings 
for — l<r<0, r = and r > respectively. The above 
discussion on kj can also be extended to all other cases. 
One can deduce that for — 1 < r < 0, < k?, while for 



r > 0, > K^. These results are numerically shown in 
Fig.(l) and the proof of these general results is the con- 
tent of Appendix C. By introducing detuned frequency, 
one might control the system ground state entanglement 
by GSI as shown in the next section. 



III. GSI V.S. ENTANGLEMENT IN DSBM: 
TWO-SPIN CASE 

Due to the fact that spins are coupled to the bo- 
son field, such interaction induces quantum correlation 
among spins. Hence, the spin system (by tracing out 



all boson states) is in general entangled. For 2 x 2 bi- 
partite system, to quantify entanglement, it has been 
proposed by Wootters0, by using the concurrence of 
the system density matrix p, which is defined by C{p) — 
Max{0, £,1 — £,2 — £.3 — £,4} where ^'s are the square root 
of the eigenvalues of spin flow matrix R defined by p: 
^ = Pi'^y ® '^y)P*{^y ® ^y)^ Subtracting in decreas- 
ing order. One should note that, after partially tracing 
out the boson degree of freedom, the spin density matrix 
belongs to the class of the generalized Werner state [T3| 
which is defined as 



PA = Trp{|j, A,/i)(j, A,/i|} 



/ Pii 



V 





P22 P23 

P32 P33 

P44 



Due to the superposition of the triplet state with j — I, 
one has P22 = P33 = P23 = P32- This density ma- 
trix has a simple formula for the concurrence C{pa) — 
2Max{p22 — v^Piipil, 0} which is determined by the com- 
petition of the populations between the entangled triplet 
and unentangled states. The form of this density matrix 
is invariant under time evolution|l4j of given in Eq.(l). 
For < K < kq, the concurrence is 



Co 



For Kx-i < K < '^A with A > 1, 



(6) 



3\/AAt^{3(4v/(l + A)3k^ + VXr^) - 4.y/3^{VTTX - VXHx - 4.ax{2VTTX - VXKl} 
9^*2(2(1 + 2A)(1 + A)k2 + Ar2) - 12^3^K^rCx - 6aA(2(2 + \)k^ - r^Xf + SafcUVSr + V^Ca) 

(7) 



where Ca = cos{^ ~ easily check that, 

when r = 0, the concurrences become 



Cx{r = 0) 



2(1 + 2A) 



(8) 



which are positive and non- vanishing for all A. There- 
fore, ground states for all k with r = are entangled. 
However, it will be shown in below that this is not true 
for r 0. Note that, for r — 0, the concurrence in- 
deed has discontinuity whenever there is a ground state 
transition(quantum phase- like transition) Q for systems 
of finite number of spins(except for N = 1). Further- 
more, in between GSI the concurrence is a constant, for 
example the concurrence for kj < k < kq is ^. The 
K-independence is due to the fact that the energy eigen- 
states is K-independent for r = 0. However, if two spins 
couple with an off-resonant mode, the characteristic of 



the concurrence is different from the on-resonant case. 
Since the eigenstates become K-dependent, the concur- 
rence is an explicit function of k. Indeed this is clearly 
shown in Fig. (2) where the concurrence of the r = 1 
case is plotted. In addition, it can be seen that for 
Kf < K < kq, one has Co > 0.5 which indicates that 
the entanglement between spins are more enhanced than 
the on-resonant result. It is also shown in Fig. (2) that 
the concurrence is strongly suppressed for Cx>i- This is 
a general tendency which also holds for r = 0. Moreover, 
as shown in Table I, Cq becomes larger as r increases. 
The suppression of Cx>i and the enhancement of Co can 
be understood by considering the ground state eigenvec- 
tors. Up to a normalization constant, states which might 
become the ground state can be expressed as follows (we 
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FIG. 1: Regions of ground states with different r. The different shadow regions correspond to different ground states. For 
example, the darkest region denotes the system ground state as |1) ect. The lowest curve is «j — ^ ii^ obtained by solving 
Eo = -f . 




FIG. 2: The concurrence with r = 1. Different shading regions correspond to different GSI. The maximum of entanglement 
appears in the interval [acj,ko]. However, for other GSI regions, the concurrence are strongly suppressed. 



neglect the labelling of j, h in the eigenkets): 
|0) = ao|0)|0)p + |l)|l)p 



2\/2k 

aA|l)|A-l)p + 6A|0)|A)p + |l)|A + l)j 



\/TTA 



(9) 



(10) 



Oo = 



|A) = 



hx = 



where A > 1. The expression of ao in Eq.(9) indicates 
the entangled triplet state has a large amplitude as r 
increases and as a result enhanced entanglement arises 
for Cq. This is due to the fact that it needs a stronger 
coupling strength to achieve GSI for large detuning(See 
Fig.(l)), therefore, stronger correlation(entanglement) 
exists. Note that this fact is also consistent with the 




results of Appendix C which requires, for large r, strong 
coupling constant for the occurrence of GSI. The eigen- 
kets for A > l(See Eq.(lO)) are different from the A = 
state(Eq.(9)) by having an extra term, namely |1)|A— l)p. 
Consequently, the existence of such term is the source of 
diluting the entanglement of the system. In passing, ob- 
serve that, in Eq.(7) the numerator vanishes by cancel- 
lation in the large A limit. It can also be understood by 
noting that, at large A, the resulting eigenket becomes: 



lA 



oo) 



{|l)| + V2|0> + |l)}|A)p. 



Therefore, the spin state inside the curly bracket is a 
separable state which implies C — Q. 

One important point should be addressed is the con- 
nection between GSI and concurrence. It has recently 
been discussed in the literatures whether concurrence is 
a good measure to quantify QPT[1|1[i1[i1[i3. For the 
case with r = 0, it has been shown that the discontinu- 
ity of concurrence is indeed associated with GSI even for 
finite system^. However, by analyzing the r = 1 case 
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TABLE I: The r-dependence of concurrence for TV = 2. 



r CTi C2_ 

-0.9 0.0977 0.0425 0.0327 

-0.5 0.3613 0.0626 0.0273 

0.5 0.0286 0.0101 
0.5 0.5691 0.0124 0.0040 

1 0.6667 0.0035 0.0008 

1.2 0.6875 0.0010 

1.3 0.6970 



"Generally, d is K-dependcnt, the data here are only showing the 
maximum values for each GSI region. This footnote also applies to 
the other tables when it is appropriate. 



TABLE IL The r-dependence of concurrence for A'^ = 3. 



r 


Co 


Ci 


Ca 


Cz 


Ci 


6 


0.5833 


0.2944 


0.0029 


0.0007 


0.0002 


7.2 


0.5942 


0.3126 


0.0017 


0.0002 





8 


0.6 


0.3233 


0.0011 








10 


0.6111 


0.3460 












carefully, one can see from Fig. (2) that there are extra 
discontinuities(See Fig. (2b)) which do not relate to GSI 
at all. They appear at the values of n where the deriva- 
tive of concurrence is discontinuous. Such non-analytical 
behavior is due to the requirement that the concurrence is 
non-negative, but not from non-analyticity of the density 
matrix. In fact this phenomenon has also been pointed 
out by Mosseri et al 15j in the model of spins interact- 
ing on a simplex embedding in a magnetic field and by 
Yangll^ for the case of XXZ chain. Furthermore, Ver- 
straete, Martin-Delgado and Cirac have also shown re- 
cently that, in gapped quantum spin system, the en- 
tanglement length is diverging without quantum phase 
transitionfivj. These results indicate that GSI and con- 
currence are not necessary in concord with one another. 
Surprisingly, a new interesting evidence is also obtained 
in DSBM. To clarify further on the relation between GSI 
and concurrence, the results of C\ for A = 0, 1 and 2 
with different r are given in Table I. From the table, it is 
noted that as r = 1.2 the concurrence of all A > 2 states 
vanishes. As r increases further, one can see more states 
possess vanishing concurrence(For example, r — 1.3 in 
Table I). Obviously, this is different from the above con- 
clusion which shows the uncorrelated non-analyticity of 
C with GSI. Here, the analyticity of concurrence is guar- 
anteed by C = 0, however, the system still shows GSI as 
K varies. In order to check this conclusion is not just the 
artifact of = 2 system, the results of A^ = 3 are re- 
ported in Table II, where pair-wise concurrences'!^ have 
been calculated. Similar to A^ = 2, for A^ = 3, one can 
see that Co and Ci are r increasing function. Moreover, 
the ground states of higher GSI with C = also appear. 
This can be seen from Table II, as r = 7.2, the con- 
currence Ci = 0, i > A. With further increasing r, more 
vanishing C; appear. This might be a finite system effect. 



TABLE III: The critical coupling constants for Ka — Kt in 
the 1-f 1 spin-boson model. 



r 








1^2 


0" 


0.7071 


0.9660 


1.4029 


1.7260 


-0.9 


0.2132 


0.2248 


0.2371 


0.2498 


-0.1 


0.4867 


0.6586 


0.9425 


1.1569 


0.1 


0.5118 


0.7043 


1.0354 


1.2758 


1 


0.5774 


0.8158 


1.2518 


1.5477 


10 


0.6770 


0.9393 


1.3910 


1.7197 


100 


0.7036 


0.9630 


1.4012 


1.7247 



" The critical {ki} of the single mode model with r — 0. 



Due to the limited spin space, after the A^th transition, 
all ground states do not change qualitatively and hence 
the ability of creating entanglement is restricted even in 
the strong coupling regime. Here, we emphasize that the 
1-1 correspondence between GSI and the discontinuities 
of concurrence seems just a special result for finite A^ 
SBM(or DM) and can not be extended in DSBM. Con- 
sequently, these results clearly establish the fact that 1-1 
correspondence between GSI and concurrence can not be 
true in general for finite A^ DSBM. It is certainly inter- 
esting to see if the 1-1 correspondence remains valid as 
A^ ^ c». 



IV. GSI WITH 1+1 MODES 

In the previous sections, ground state instabilities of 
DSBM and its correlation with concurrence have been 
treated. It is interesting to see if the conclusions still 
hold for the case of multi-mode spin-boson model. In this 
section the case of two-mode (1+1 mode) model will be 
analyze by adding one off-resonant mode to the resonant 
SBM. The Hamiltonian is 

H = Ha + Hi 
Ho = Jz + a^a + b^'b 
Hi = rb'^b + HaJ+a + nl^J^a^ 
+HbJ+b + nlJ^b^ 

where {fo, 6^} are annihilation and creation boson opera- 
tor of the off-resonant mode with ki, being the coupling 
constant. Similar to the single mode model, it is easy to 
check [Ho , Hj] — and one can choose the eigenstates of 
Hq to represent operator Hj: 

Ho\j,m)\na,nb)p = X\j,rn)\na,nb)p 
Hi\\h) = h\X,h) 
H\X,h) = Exh\X,h) 

where A = rta + rtf, + to and Exh = X + h. The dimension 
of Hi is 1 for A -1 and 3(A + 1) for A > 0. Similarly, 
we will omit the labelling, j and h, of the eigenkets in the 
following discussions. GSI of the 1+1 model with r = 1 
is shown in Figs. (3). For example, by fixing Kb = 0.4 the 




(c) 



Kb ^ 




-1 




1 1 S 9 






FIG. 3; The eigenenergy of A = -1, 0, 1, 2, 3 for r = 1. (a) = 0.4; (b) Ka 
in (c) and (d) indicate the eigenstate with A = and A = 1 respectively. 



0.4; (c) Kb ~ 1.1; (d) Ki, = 1.8. The black arrows 



sequential GSI, namely {Ei ^ i?o — *■ i?i — *■ i?2 ^ • ■ •} 
are clearly shown in Fig. (3a) ( The result of keeping 
fixed at 0.4 is shown in Fig. (3b)). However, as Kb in- 
creases, the state |1) = |T)|0a,0f,)p may no longer be the 
ground state. Therefore the sequence of GSI does not 
have to begin from A = — 1. This result can be seen in 
Fig. (3c) which shows the sequence of GSI with ki, — 1.1 
as {Eo El ^ E2 ^ ■ ■ ■}■ This fact just reflects the 
result of GSI in DSBM discussed in Sec. III. Moreover, 
if Kb is further increased, more lower spin sector eigen- 
states get kicked out of the GSI sequence, this is shown 
in Fig.(3d) with Kb = 1.8. 

It is also interesting to see how the detuning parameter 
is related to the pattern of GSI. By adding one extra off- 
resonant mode to SBM, the results are shown in Table 
III which also includes the results of SBM for compar- 
ison. To keep things simple, the results are evaluated 
with Ka = Kb without losing generality. One can see that 
the critical values (k^) for GSI are increasing function of 
r but being bounded by the results of SBM(See the first 
line of the table). Similarly, for r closed to —1, Ki also 
approach the results of SBM which can be determined 
explicitly from Fig.(l). These interesting results can also 
be understood from the energy spectrum. For example, 
consider kj which indicates the crossing of the eigenen- 
ergies of |1) and |0). When r is large, the off-resonant 
mode has higher energy and its excitation costs more 
energy. Thus, |0) involves dominantly the lower energy 
boson which is the on-resonant mode. As a consequence, 
the determination of kj is governed by SBM. Physically, 
what is happening is the effect that is well-known in most 
physical systems, namely, the result of decoupling effect 
of far off-resonant driving. On the other hand, at the 
limit of r — > — 1, the off-resonant mode with lower fre- 
quency is dominating. Therefore, |0) can have more con- 
tributions from the off-resonant photon and consequently 
the determination of kj is dictated by DSBM. It is impor- 
tant to point out that with extra mode, Ki can be reduced 
significantly. For example, as r = 1, = 0.5774 which 
is smaller than 0.7071 and 1 of the critical couplings ki 
with r = and r = 1 respectively(See Fig.(l)). Since for 
K > Ki the ground state is entangled, it is certainly im- 
portant to obtain GSI in the weak coupling regime. How- 
ever, having GSI at lower critical coupling is not enough 
for practical reasons. One important requirement for em- 
ploying entanglement in quantum information science is 
to have strong enough entanglement or maximally entan- 



TABLE IV; The r-dependence of concurrence in 1- 
boson model. 



Co 



1 spin- 



Co 





-0.9 
-0.1 
1 

1.8 

5 
10 



0.5 
0.1074 
0.4898 
0.5455 
0.5462 
0.5381 
0.5253 



0.0286 
0.0508 
0.0324 
0.0211 
0.0226 
0.0316 
0.0348 



0.0101 
0.0394 
0.0114 
0.0075 
0.0084 
0.0135 
0.0161 



gled state. Hence, it is necessary to see if adding extra 
mode can either enhance or suppress entanglement. 

The r-dependent results of few C^'s are tabulated in 
Table IV which also contains the results of SBM (the first 
row). For practical aspect, we only concentrate on Cq 
which has higher entanglement. One can see that the 
maximum value is around r — 1.8. However, at larger 
r, Co is decreasing toward the value of the SBM result. 
This result is the same decoupling effect discussed pre- 
viously and once again the on-resonant mode determines 
Cq. Note that the maximal value of Co is higher than 
the one in SBM and the determination of the value of 
r with the maximal concurrence is a balanced result of 
the competition between two modes. The same effect 
happens for other Ci in Table IV. One should note that 
the entan glem ent obtained in the weak coupling region is 
distillable[l9j and can be enhanced comparing with the 
on-resonant result(C = 0.5). Therefore, one can have an 
" entanglement switch" Q by controlling the first ground 
state transition and it seems that having extra mode can 
do just that. All in all, this result seems to suggest that 
1-1-1 mode might be functioning better than mono-mode 
models. It is important to justify the above conjecture 
by studying a system with three cavity modes which will 
be reported elsewhere. 

Finally, for comparison, the results of both cavity 
modes being off- resonant with ra = 1.2 and rf, are 
given in Table V. The result is quite interesting since 
the value of Co can be higher than the corresponding 
results(r = 1.2 with Co = 0.6875) of SBM in Table I. For 
example, when rb — 1.5 one has Cq — 0.6998 which is 
larger than 0.6875 given in Table I. Furthermore, as rf, 
becomes very large, one can see from this table that Co 
approach the value in SBM. Again this is just the effect of 
decoupling. Further support of this fact is revealed the 
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TABLE V: The results of 1+1 cavity modes with ra = 1.2 

and Ka = Kb- 



n 


Co 


Ci 


C2 


1 


0.6764 


0.0023 


0.0004 


1.1 


0.6823 


0.0016 


0.0002 


1.3 


0.6921 


0.0005 





1.5 


0.6998 








100 


0.6920 


0.0043 


0.0035 


10000 


0.6875 


0.0011 






values of Ci and C2 in Table V. Hence, for enhancing 
entanglement, it is likely that off- resonant multi-mode 
model is a better candidate and deserve further analysis. 
Similarly, the absence of discontinuity of the concurrence 
with GSI can also be found in the multi-mode system (Sec 
Table V). Moreover, we find that the extra discontinuities 
of concurrence also appear in two off-resonant modes of 
spin-boson model. These facts show more supports for 
disconnecting GSI and discontinuities of concurrence in 
spin-boson model. 

V. CONCLUSION 



a continuous concurrence. Even though the above results 
are obtained for N = 2 system, we have also obtained 
numerical results for iV = 3 which also support our con- 
clusion. Furthermore, the 1+1 mode model is analyzed 
and interesting results are obtained. It is seen that the 
effects of extra mode are two folds. First of all, GSI can 
happen at weak critical couplings which are important 
for having entangled ground state. Secondly, the entan- 
gled ground state with extra off-resonant mode have a 
higher concurrence comparing to the results of SBM and 
DSBM. These results may be useful in the context of 
quantum information science. Moreover, we also con- 
firmed that GSI and the discontinuities of concurrence 
remain uncorrelated even extra mode is included. There 
are several directions for further study along this work. 
It is interesting to obtain results for adding more modes 
to justify the effects of extra modes as obtained here. 
It is also necessary to analyzed the A'' atoms case. Fur- 
thermore, results for more than two-level system are very 
important. For example, for 3-levcl systems, it is inter- 
esting to analyze either the A system or V system to see 
if the results obtained in this work remain valid, since 
these systems arc also quite common in atomic physics. 
These problems will be pursuit in the future. 



In this work, we show explicitly that GSI and entan- 
glement are not necessarily connected. This is shown by 
a rigorous analysis of the spin-boson model(SBM) with 
two spins. By knowing the full spectrum, it is shown 
that the sequential quantum phase transitions occur in 
this system and the closed form expression of concur- 
rence is obtained. Employing these results we are able to 
clarify the the relation between GSI and entanglement. 
Contrary to the results in literatures which are concen- 
trated on spin-spin interating systems, it is shown clearly 
that in the detuned spin-boson model, concurrence is not 
a good measure for quantifying GSI. This is shown by re- 
alizing that not all the discontinuity of concurrence are 
associated with the ones appearing in GSI and on the 
contrary the system having GSI can be corresponding to 



This work was supported by the National Science 
Council of R.O.C. under the Grant No. NSC 93-2112- 
M-006-006 and 94-2112-M-006-016. The support from 
the Center for QIS at NCKU and the Tainan Branch of 
the National Center for Theoretical Sciences of R.O.C. 
are acknowledged. 



APPENDIX A 

In this appendix, we provide the finite matrices of Hj 
for given A with j = For A < the general ma- 
trix form of Hi for arbitrary A^-atom system is: (with 
excitation number X = + u and v < N) 



\ 







kr K*Vk + \ ^Rv-k 

K^fl7^^fB^ {y - l)r K*-sfv^Rl 

V K^\fRx vr ) 

(Al) 

where = x{N + 1 .x) and i?2;<o = 0. The non-vanishing off-diagonal elements only exist next to the diagonal on 
each row. For A > the dimension of the matrix is fixed which is {N + 1) x [N + 1). The general matrix form of 



/ k*s/r;^ 

K^/R^ r K*^/2^R^^i 

kV2^/R~[ 2r K*V3y/R^ 



K\fksJ Rv-k-\ 
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Hi is 



V 



(A2) 



where \ = ^ +v and i/ > 0. For the first ground state 
transition, the critical coupling for arbitrary A'' and de- 
tuning r is 



1 + r 
N 



(A3) 



which can be obtained by solving the equation of Es = 



En_,-, where 



En = — - 



AT 



(A4a) 



En^i = -y + l + ^{r- V4iWH^}-(A4b) 



APPENDIX B 

This appendix is to prove that not only the existence 
of GSI of two-atom system with arbitrary detuning but 
also these GSI occur in sequence. In order to prove these 
results, it is necessary to determine the ranges of (fix and 
6* A- Denoting x = (j)"^ and recalling the definition of (p\, 
one has: 

<px = cos { ^=^1 



'at 



cos ^{ 



sen (r)3\/3a; 



} (Bl) 



where ta = 1 + 2A and sgn(r) is the sign function. Since 
< K < oo and — 1 < r < cjo, one has < x < oo. The 
ranges of (fix and 6x are: 



cos ^( — ) < V'A < cos ^( — -) 
rx Tx 

Ox<Ox< Ot 



-{tt-cos ( — )}. 

O TX 



(B2) 
(B3) 

(B4) 



Although the maximal value of 9x is A dependent, it is 
easy to check that 6x is bounded as follows: 



(B5) 



For (i), it is easy to see that Eq is a monotonic de- 
creasing function by directly checking 9„£'o: 



d.Eo = 



4k 



(B6) 



which is negative definite for k > and — 1 < r < oo. For 
Ex given by Eq.(5c) with A > 1, the partial K-derivative 
of Ex is: 



d.Ex = 

nx = 



(B7) 



TX 



3ax\/a] 



I - 27K4r2 -I- ^/3r{TxK^ - r^) smOx. 



Note that, a\ > TJk^'^ and ax > 0, then f^A > for 
all K and r. Moreover, by using the ranges of the angles, 
one has 

0.804 ~ cos 6*+ < cose'A < cos6'i' ~ 0.917 
0.399 ~ sin6'i" < sin6'A < sin 61+ - 0.595 

Therefore, cos^a > sin^A > 0. For r > and ta > 3, 
one has Xa > Xa - Xb with 



XA = \ 3{al-27K'^r^)cos^ Ox + 3K^r sin t 



XB = \l^r^ sin^ ^A, 



By expanding out ax in xa and regrouping terms one 
obtains 

Xa = \J 3r^ cos^ 0a + ^ + 3K^r sin Ox 

where A denotes the remaining positive part inside the 
square root. Obviously, cos^ 6'a > svc? Qx, one concludes 
Xa > Xb =^ Xa > for all k and r > 0. Similarly, it is 
easy to check, for — 1 < r < 0, xa > for all k. Therefore, 
Ex are monotonic decreasing function for A > 0. 

To show f{K,r,\) = -Ea+1 ^ Ex is a monotonic de- 
creasing function of k for any r and A, it is necessary to 
break down the proof for different regions of A. For the 
A = — 1 case, by Eq.(5), 



dJ{K,r,l) = d,{Eo-Ei) 
= d^Eo 



(B8) 



By knowing the ranges of the angles, we are now in the 
position of showing the conditions (i), (ii) and (Hi). 



which is monotonic decreasing as proved in the criterion 
(i). One can easily check /(k, r, 1) approaches 1 at small 
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K, while becomes — oo at larger k. Therefore, the critical 
coupling ki determined by /(k, r, I) = uniquely exists. 
For A > 1, one should show that the k derivative of / does 
not change sign for all k. Alternatively, it is equivalent to 
show the function g{K, r, A) = d^Ex being a monotonic 
decreasing function in A, such that it ensures g{K,r,X + 
1) - gin, r. A) = d^Ex+i - d^Ex = d^f < 0. By using 
Eq.(B7), 



a.,(A, .,r) = 3x3.Ex = -^^rj=^^TW^ 

{V3{SK^r + ^Jy^) COS0X - 3(K^r - ^Jy^) sIu^a} 



X 



where 



3(16A2 
2(16A2 



16A 
16A 



For r > 0, k r < yni, one has dcg{X,K,r) < 0. Simi- 
larly, dcg{X, K,r) < is still true for — 1 < r < due to 
3K^[r| < s/y^. Thus one has dcg < for all k, A and 
r. Therefore, we have shown g{K,r,X) is a monotonic 
decreasing function of A and then /(«;, r, A) is a strictly 
decreasing function of k. Furthermore, one can check 
that, for A > 1, 

/(K,r,A)U^o = l + r>0 
/(K,r,A)U_>oo = V2T4A-V6T4A<0. 

As a result, the crossings {k,i}{i > 1) has unique solution. 
To prove the remaining case with A = 0, we express 
r, 0) in terms of x: 

dJ{K,r,0) = d,{Ei-Eo)=vT 

4 

V = 



V'(8a;+l)(6a; + l) 



F = y/xi6x + 1) - y/3x{8x + 1) cos 6*1 



-sgn(r)(3a; - !)< 



Ix{8x + 1) 



sin 6i 



where y = y\/r^ . If F is negative for all x > 0, then 
/(k, r, 0) is monotonic decreasing. Let us start with 
sgn(r) = +. For x > \, V6x + 1 < V&F+T and 
•\/3cos^i > 1, then ^/x{Qx'+Y) < x{8x + l)-\/3cos0i. 
Therefore, we have F < for a; > |and r > 0. Similarly 
one has F < for .t G [(0, \) with sgn(r) = — . However, 
for X e [0, ^] with sgn(r) = +, the last term of F is neg- 
ative, thus it is not obvious that F is negative definite. 
Therefore a different approach is called for. To proceed 
further for x G [0, i] with sgn(r) = -|-, let 

A = yx(6^TT) 



„ , xi8x + 1) . ^ 
C = ^/3x{8xTTjcos6l, 



and then T = A + B ~ C. Our logic to prove F is still 
negative is to show that there exists a. 6 > such that 
SC>A and (1 - 6)C > B, then C > A + B <Oioi 
X e (0, 5). To begin with, one has 



50A 



The other condition is 




6x + l 
3(8xTT) 



sec 61 



(B9) 



5 < 1- 



1 



3x ^ 
— tan 61 



/3y 



(BIO) 



Combining Eq.(B9) and (BIO), 
1 



1 - 



3a; 

— tan fill >6> 



6x + 1 



sec 01 >0 



/3y ' ~ ' ~ V 3(8-^ + 1) 

Therefore, 5 exists if the following condition is satisfied: 



,1-3^, „ 
1 — I — ■^==- tan 61 



/3y 



6x + l 
3{8xTl) 



sec6'i}>0 (Bll) 



Z2 



It is easy to show the ranges of Zi and Z2 are: 



0< <- 



— sec 6*5^ 



Therefore, 
2 



<Zi + Z2< - 

3 - 3 



sec 9+ < 1 



(B12) 



By the same approach, one can show that F is negative 
definite for r > | with sgn(r) = — . Furthermore, it is 
easy to show that /(K,r, 0) 1 when k is small and 
changes sign at large n. This completes the proof of 
showing K2 exists. For now, we have shown the crossings 
{Hi} between Ex and Ex+i exist even in the detuning 
two-atom system. In what following, we will show these 
crossings occur in sequence. 

We start from A = — 1 case, with the solution of Ei = 

Eq, Hi = \/^-^, one has: 



^1 - ^0 = (2 + r) - -^V(2 + r-)2-(l + r)cos^i 
v3 

where the tilde symbol denotes quantity at the appropri- 
ate critical k, for here it is k = kj. For A = 0, imposing 
El = Eo, one obtains 



1 + r 
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And 



E2 — El 

2 / 

1 + r - -{v3q;2Cos^2 
o 



'3ai cos ^1} 



r 4-\/ 6ko^ + r"^ cos^i , 1 / IOkq + cos 6*2 
2^ ^ ~ 2Y 6K2 + r2 cos^i 



\/3 Skl + r'^ 1 



(B14) 



8 Y 6Kg + r2cos6ii 
where Eq.(B13) has been used in the second line and 

2 A, _i 3V3Kl+r) 

cost*! = COSj-ITT — cos )[ 

n A I -1 3\/3r(l + r) 

COS6'2 = COSj-lTT — cos ) \. 

^3^ 2V(5 + 5r + r2)3^^ 

In order to estimate the value of Eq.(B14), By using the 
inequalities (B3) and (B5) one has: 



cos 



^0+ 



cos 6^ 
1 



< 



COS6'2 



^ < 



COS fi 



< 



COS ^2 

COS 9f 

1 



COS^i 



For sgn(r) = +, it is easy to determine the minimal value 
of the part in the bracket: 



{•••} > 1 



1 /10k2 +r2 cos6'^ 
2y 6kl + r^ cos 61+ 



V3 8kl + r^ 1 



8 y GRl + r^ cos 0t 



V3 /I 1 



2 V 3cos6»5^ 8 V 3cos6»^ 



(B15) 



Then, E2 > Ei at Kq for r > 0. For A > 1 cases, from 
E\ = Ex+i, 



1 



Therefore, 



-{^3aA+icos^A+i - ^/3de^cos0x}- (B16) 



{Ex+2 - Ex+i}\ii, 



= {r + 1) - -{^/Mx+2Cos6x+2- \^Mx+icos6x+i} 



= -{2a/3q!a+i cos^A+i - \/3aAC0s6'A 
-a/35^cos^a+2} 



(B17) 



where Eq.(B14) has been used in the second line. For the 
ease of discussion, wc denote ^(A) = \/3ax cos Ox- If the 
curvature of ^(A) is negative, then we have 2^(A + 1) > 
^(A)+C(A+2) which is just the condition ofEq.(B15)> 0. 
By taking second derivative with A directly, we obtain 

3WA) = -(2i)i(^)x 

{V%acos6'a + 27^ sin6iA} < 0. (B18) 

Therefore, the crossing {/tA}|A>i actually occur in se- 
quence for all — 1 < r < 00. This complete the proof for 
having sequential GSI. ^ 



APPENDIX C 

In order to prove Kf < < Kf where nj^'^'''^ are for 
r{<, —, >}0 respectively, the hamiltonian of the system 
can be rearranged as: 



H° + ra^a 



(CI) 



where denotes the on-resonance hamiltonian and 

is the detuned hamiltonian. a^a is photon num- 
ber operator which is positive-valued. For r > 0, given 
K = nf, the ground state of the system is Ig^^^i) such 
that H'^\gf'j^_i) = E^\g^_^_^). The expectation value of 
in Idi+i) '^^^ written as: 



E^ = {gr+,\H^\gr+,) + r{a^a),> 
> {g'W)U,>+r{a^a)g> 
E'>E''\.. 



(C2) 



Where |.g°) is the ground state of H'^ . The inequality 
occur by noting that \gf+i) is a trial state for . More- 
over, due to the fact of sequential GSI proven in appendix 
B, Eq.(C2) implies k° < nf . Similarly, for -1 < r < 0, 
instead of using \g^j^i), we evaluate the expectation value 
of with the ground state of at k = and 

obtain the desired result : 

E"" = (ff°+i|i?l5°+i> + |r|(aM,o 
> E'-\,^,o + \r\{a)a)gO 

£;n«=«o < E°. (C3) 

Again the property of trial state has been used in Eq.(C3) 
which implies < k^. ^ 
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